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Abstract 

Charge and spin-orderings in the i-filled organic CT soUds are of strong interest, especially in view of their possible 
relations to organic superconductivity. We show that the charge order (CO) in both ID and 2D CT solids is of 
the ...1100... type, in contradiction to mean field prediction of ...1010... CO. We present detailed computations 
for metal-insulator and magnetic insulator-insulator transitions in the S-ET materials. Complete agreement with 
experiments in several 9 systems is found. Similar comparisons between theory and experiments in TCNQ, TMTTF, 
TMTSF, and ET materials prove the ubiquity of this phenomenon. 
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Ground states involving spatially inhomoge- 
neous charge distributions have in recent years 
been observed in a wide range of novel electronic 
solids. In the case of the ^-filled band organic 
charge-transfer solids (CTS), it has been suggested 
that the appearance of charge-order (CO) may be 
related to superconductivity (SC) [1,2]. Two pos- 
sible CO patterns have been suggested in ID sys- 
tems, corresponding to "...1010..." and "...1100..." 
site occupancies (see Figs. 1(a) and (b)). Here '1' 
and '0' correspond to actual electron or hole den- 
sities 0.5-l-e and 0.5-e, respectively. The ...1010... 
CO corresponds to a 4kF CDW, where = 7r/2a 
is the Fermi wavevector for non-interacting elec- 
trons. In the presence of strong nearest-neighbor 
(n-n) Coulomb repulsion V, the obvious preferred 
state is ...1010... [3]. However, this CO can occur 
only above a V > Vcr (see below). Furthermore, 
it is not clear that this CO can explain the full 
T-dependent behavior of the ^-filled band CTS. 
The ...1100... CO pattern occurs for V < Vcr, and 
is a consequence of cooperative electron-electron 
(e-e) and electron-phonon (e-ph) interactions [1]. 
This cooperation may be the key element behind 
the unconventional SC in the CTS. 
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Fig. 1. (a) ...1100... CO state. Filled (empty) circles corre- 
spond to sites with more (less) charge. Lines indicate bond 
order strength, (b) ...1010... CO state below spin-Peierls 
transition (c) The BCSDW in weakly 2D. 

Here we examine the theoretical and experimen- 
tal evidence for both ...1010... and ...1100... CO 
in the quasi- ID, weakly 2D, and isotropic 2D lat- 
tices. We conclude that with the possible excep- 
tion of (TMTTF)2X, in all cases the pattern of the 
CO is ...1100... Further theoretical and experimen- 
tal work is necessary to understand (TMTTF)2X 
completely. A second important issue we examine 
is the relative temperatures at which CO, metal- 
insulator (MI), and magnetic transitions occur. 

In ID, we consider the Peierls extended Hubbard 
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Fig. 2. (a) Schematic structure of the (?-(ET)2X conduct- 
ing layers, consisting of the ET cationic molecules. Solid 
lines correspond to stronger hopping tp, dashed lines to 
weaker hopping tc- (b) vertical stripe, (c) diagonal stripe, 
(d) horizontal stripe. Filled (unfilled) circles correspond to 
molecular sites with greater (smaller) charge density. 

Hamiltonian, H = Hq + Hee , where 
Ho =- ^[io - a{uj+i - Uj)][c]^Cj+i + h.c] 

jo- 
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In TJq J a and /3 inter- and on-site e-ph couplings and 
K\, K2 the corresponding spring constants. The 
displacement of molecular sites from equilibrium is 
Uj , and the on-site phonon coordinate is Vj . U and 
V are on-site and n-n Coulomb repulsions. 

It is well known that for a = /? = 0, the ...1010... 
state occurs only for V larger than some criti- 
cal value, Vcr- For V < Vcr, the Luttinger liquid 
state occurs, which becomes the ...1100... state for 
nonzero a or /3[1]. In the limit U —>■ 00, Vcr — '2to, 
and for finite U, Vcr > 2to [4,5]. Consequently, 
for realistic U, only for systems with the smallest 
to there exists a narrow range of V for which the 
...1010... is a real possibility [6]. Since precise es- 
timation of V is difficult, an alternate approach 
to determine the CO pattern is to examine the 
bond distortion pattern below the spin-Peierls 
(SP) transition [6]. A SP transition within the 
...1100... CO has bond pattern SWSW (S=strong, 
W=weak,W' t^W; see Fig. 1(a)), while the bond 
pattern for the SP phase with ...1010... CO is 
SSWW [6]. Experimental observations of SWSW 
in MEM(TCNQ)2 and other ID 1:2 TCNQ sys- 
tems [7,8] provide direct evidence for ...1100... in 
the 1:2 anionic CTS. In (TMTTF)2X with smaller 
to (larger V/to) CO appears below the MI transi- 



tion and the situation is more ambiguous. There 

exists currently no information on the bond dis- 
tortion pattern in the SP phase of (TMTTF)2X. 
It has been suggested that there exists a phase 
boundary between the CO and SP phases in 
(TMTTF)2X [9]. This might indicate a switching 
from ...1010... at high T to ...1100... at low T [10]. 

In the (TMTSF)2X materials, CO coexists with 
SDW order of the same (2ki?) periodicity [11], 
which indicates the coexisting bond-charge-spin 
density wave (BCSDW) state (Fig. 1(c)) that oc- 
curs within a weakly 2D model consisting of ID 
chains with weak interchain interactions t_L [1]. 
Evidence for the BCSDW (with very weak mo- 
ment) is also seen in a-(ET)2X [1]. 

We next consider 2D 6'-(ET) lattices (see Fig. 2). 
CO here corresponds to "stripe" patterns shown in 
Fig. 2 [12]. In the Hmit of a = /3 = the 2D Hamil- 
tonian includes weak and strong hoppings tc and tp 
(see Fig. 2(a)) as well as Coulomb interactions C/, 
Vc, and Vp. For realistic systems Vc ^Vp = V, and 
the classical energies of all three stripes in Fig. 2 
are equal. Within Hartree-Fock (HF) theory [12], 
the vertical stripe state is found to dominate. 

We have performed exact diagonalization stud- 
ies for a 16-site 0-latticc with 8 electrons, periodic 
in all directions (see Fig. 4). Instead of explicitly 
including the e-ph interactions as in ID, we follow 
our previously established procedure[l] and add a 
site energy component ^ - ani , where is negative 
(positive) for the '1' ('0') sites. This is equivalent to 
including /3, with fixed Vj . In Fig. 3 we plot the en- 
ergy gained upon stripe formation, AE = E{ei = 
0) — E{ei 7^ 0) for the all three stripes against V , 
for U = 0.7 eV, tc = -0.01 eV, and tp = 0.14 eV, 
for |ei| = 0.01 eV. We find that the vertical stripe 
is never the ground state and the horizontal stripe 
dominates for all V > 0.18 eV. Importantly, while 
the vertical and diagonal stripes are ...1010... along 
two of the three directions, the horizontal stripe is 
...1100... along the both p directions (see Fig. 2). 

We also calculate the charge densities rij and 
bond orders, = X^^((c|'^Cja + cJ^q^)). The 
bond orders are spontaneously distorted for the 
horizontal stripe state in spite of uniform tc and tp. 
This implies spontaneous c-direction dimerization 
and p-dircction tctramerization for any a = 0+ 
(see Fig. 4). No such bond distortions arc found in 
the vertical or diagonal stripes. 

The above ground state calculations arc rele- 
vant only at low temperature. Experimentally, CO 
at high T (~190 K) [13], and a spin gap at lower 
T (< 50K) [14] have been observed. Transitions 
at finite temperature are determined not by the 
ground state energy, but by the free energy. At tem- 
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Fig. 3. The energy gained upon stripe formation, for the 
three stripe patterns of Fig. 2. For V > 0.18 eV, the ground 
state has the horizontal stripe CO. 



peratures T < Tmi , the excitations of a strongly 
correlated system are predominantly spin excita- 
tions. Because of the greater multiplicities of high 
spin states, the free energy at high T (but below 
Tm/) is dominated by high spin states. We there- 
fore perform similar calculations of nj and bij for 
the ferromagnetic (FM) spin configuration in the 
horizontal stripe (total spin S = Smax = 4). The 
bond dimerization in the c-direction persists even 
here, although the p-bonds are dimerized and not 
tetramerized. The interpretation of the S = and 
the S = Smax calculations, taken together, is as 
follows. As T is lowered, formation of the horizon- 
tal stripe CO occurs as soon as charge excitations 
are too high in energy and the free energy is domi- 
nated by high spin states. The CO is accompanied 
by c-axis bond dimerization, in complete agree- 
ment with experiments [14]. At still lower T, high 
S states are thermally inaccessible and the free en- 
ergy is dominated by the S = state, where there 
is a second transition, also in agreement with ex- 
periments. We believe that the second transition 
involves tetramerization of the p-bonds. We have 
argued elsewhere [15] that this change in the bond 
distortion pattern from dimerization to tetramer- 
ization is accompanied by a spin gap[15]. We note 
that from analysis of ^^C NMR lineshapes Chiba 
et. al. have concluded that the CO in 6-(ET) cor- 
responds to the horizontal stripe [16]. 

In conclusion, the dominant CO pattern in the 
i-filled band CTS appears to be ...1100... . Direct 
evidence for this exists in the 1:2 TCNQ (SWSW 
bonds below the SP transition), (TMTSF)2X 
(mixed CDW-SDW with same periodicities for 
CDW and SDW), and 6l-(ET)2X (c-axis bond 
dimerization at Tco and the occurrence of two dis- 
tinct transitions). Further experimental work, in 
particular elucidation of the bond distortion pat- 




Fig. 4. The periodic 16 site lattice investigated for all 
three stripe CO orders. Filled (open) circles indicate large 
(small) charge density and numbers are bond orders. See 
text for parameters. Note the bond dimerization along 
the c-direction and the bond tetramerization along the 
p-directions. 

tern below Tsp would be necessary to determine 
the CO in (TMTTF)2X. Most importantly, the 
occurrence of the ...1100... CO is a signature of co- 
operative e-e and e-ph interactions and this raises 
the interesting possibility that the unconventional 
SC in CTS is a consequence of this cooperation. 
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